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ABSTRACT 


Algorithms for resolving the point-in-polygon problem 
are. surveyed. The scope of the survey includes only those 
algorithms, and consequently those applications where 
polygonal vertices are represented using vector rather than 


raster notation. 


Two classes of algorithms are described. TMhessdirect 
methods sequentially process the n vertices of the polygon 
and require O(n) operations to resolve the point-in-polygon 
problem. The pre-conditioned methods operate on a 
hierarchical preprocessed representation of the polygon, and 
anemwOrmcomplexity O(log nn); O(log2n), Or O(n). sHOowever, ecoe 
pre-conditioned algorithm have associated with them 


additional preprocessing and storage costs. 


iv 


n oy : 
Tanned: 

pany 
: — 
¢ 


nat? a$afex rotvev Shiau as ea ik Sanrsiey | 


jrebaese ede F 


_ 


ifeslt, oat 4 ghodisass6 stn saisixeprs Zo 2ueasl> ObF . 
nogvdey ot. to esuidvev.o ‘ent: apsdes; ylisisteupes eedaee | 
Adwv ien=ti-$iing-.922 ‘eicess Stato se qeno (nd0 etiepe2 bas 
6 me 363500 eberitan, (bencigiands-osq) - )Sn? melon - ; 
bas ,rdeyiog Sits «2p 10 LegRe2O 1023 Poteser westg ievitozesea@) 
07 .govarcH «ino: ai taser io" mat elie yrSk 3 Lures 20 8 


ora 12 41.) is tet4oces aye ig soplS bent it Lbeoees 


p39) soexotd Bie len bese oigeig ienoldibte 7 
iY 


ACKNOWLEDGEMENTS 


TOmMVEESUDGCLV 1 SOL, Dig tocanmcCabay, Sincere thanks for 
the patience and guidance he provided me during’ the 
research. The thesis committee members Drs. Wayne Davis, 
Jean Claude Muller and Len Schubert read and commented on all 
chapters of this thesis. 

I must mention my good friend, Kofi Marfo whose 
invaluable help enabled the generation of numerous diagrams. 

My sponsor, the Department of Computing Science, 
provided the essential financial assistance. 

Dr. Narh Omaboe of Trenton, Ontario, gave all the moral 
SUDDOREMMeeOeC ll BEL On OLIN es aDOULMmLO LSE meCOMDLOtemmmIeSCarCc) 
document. 


My sincere thanks to them all. 


: 


107 aha ear “wade ‘test ad eit le yor 7 . n 
ere aay 


is 
4 


ime 


sit pais a0 ina Selig ay sonabiug tos 


ap 2 : 2 -outiatenm Bleeds gat — 
,eivet sovr¥ «270 aretmem pase Li é 


—o . 
shun iD wast 


— 

= 

or 
om 
- 


Lis no Setecwntes tes hes 7 fudat oo! 


ateoda aid? 30 woetqeia 7 


seoilw oftewet Flos ,b608T2 SA0E Yr san saaum 1 : o 
eeewk bib anctsdan 20 mebiSsencs otis (iceste GI 4 etdsviaval 7 
. ayAes Se on! Pum ~* Sheen! Tce. 0 ‘7 i os To pe va 
Sibdassee | ns cl 23 ads tenivesq 
- 
attyi oa ie _ aE - =~ ine? fre 2 Aw’ 4 anh @234 ,3a 
- ~~ 
MISHSe as a5 tm win 105% sniad . o2 yaheen + cope” 


/2eecw2ob 


Chapter 
DNR DUC LON @me teeters erategsier siete 


Is 


ie 


IV 


al: 
12 
as 


Table of Contents 


BaCRORO UNC acre: crete le eres ts 


Statement of the Problem . 


Outs Tne FOEwTh eC Saks) eae. 


DERE CilmME TODO Merete steleistelecre 


2 ol 
22 
2e3 
2.4 
2a 
2.6 


iBpacclele(bletesmele’ 4 Ser 
Directed Line Method 


Area of the Polygon Method 


Sum of Angles Method 
Line Crossing Method 


Summary eoeeeeeege#eee#ee#eee 


PRE=CONDITIONED METHODS .. 


Slealt 
She 
5 f:8) 
3.4 
3.) 
3.6 
Brea! 


Mr EieO cd U Cia1 Olumerersucrsrer sss 
Sieueigh Wika eleer wg Aa Aaa 
Chatn-Method).... 1... aaete 
Talia UCmMeC COG. artes 
Rectangle Method .... 
Seenke, Welles so doen oo 


CONC WITS ONaesters sretererel ete 


PRACTICAL CONSIDERATIONS . 


4.1 
4.2 
4.3 
4.4 


CONCLUDING REMARKS 


Introduction eoeee#e#ee#ee? 


20: 0. 8.8) 


eee37e#e 


Rectangular Enclosure Method ... 


Circular Enclosure Method 


eoeee#ee 


Multiple Polygon Processing .... 


eeee0no5xoe3eeeeeeete@e¢@ @ 


(eo) Sp en fey) 


He 
14 
ie) 


Zu 
Zu 
22 
28 
40 
oul 
64 
82 


83 
83 
83 
84 
86 


88 


ve 
8 ndode mw eee eees Ben ee orheit again * | 
6 a Pee Me P| ae 4-ues it ssiaiers Jak L 
if Oe Aut tte. afi fil. ‘bedso7a a. 
2 seO enn enka se aiaine fSie ae’. meas : ect2. 20 Sook , 


ti aeibprenenses eee e Hopasn, aha to Ge res ; 
>i «gp bt eee de Seed enet= oe Rowen on kerebaS sqra é.t 
: +) yang ae 


- FL sab aon apes dS Ne tne me she a babes 


ig 988 a oo. 04 we ms =~ 204 eee SOUL SSH CSD Res A | 
is A a Se SE elas AVL say oratn >. ree 
ss enc GS teh Pate eee nae tap <1 es <r forltg4 PAAW> a 


CF “ake stn aber ar se as ert« cheeks’ « foraem Arete -24E 
en rr er ee ee Boaseh olpaagat ‘bee 
TE ew dd doa br el ee ereoeratyeegs DUG IIE a kl tt 
ba Geieshidse PAA Shdes a 1. 2Soreaee Ge 
ou OS te Kae oe a id Rea ho et=as) Vv.e 


— a) 
7 ‘ 


ea sean bed een ae * eee oper EiOLfansazeue? JAS L?(3an9 
- as Samed tee ots eee Enos ee er eee ho2rgubled tat ink 
Deir -pitidaree fodteh euorotoct telugmadoet S23 
ant se aaessyhctes =; forsan Pulleeiogs 26.05 255 i.e : 
Maen peastren >> Seen noygio® SigkatuN th © 
a 
Dae ite ee 
ne eich gas bpraciaae BAI sate hd . 


> 


LIST OF TABLES 


Table Page 
I SHMMENEN, Cone IOivacvels Mines Ania ro on Gono OD do AS OS iy) 
Tt Experimental Results using a Polygon of 

Oy OMe Cia © SuMsteretel ets chelsisteroteieretslcrscete averse, ste) etait tsts ets 1 
Teel SUMMaGy POLE Lem CONG Le 1 ONCGUMMe LNOGUSa mts clniel siete le 82 


aed 


tn = iw fi Legis Fangalt 
sane eve sei apres . aah — 


; ~ 


LIST OF FIGURES 


Figure Page 
Teo Positions O0 “a "given Point 2=(x,y) 

Rea GiVviegrO Fae? Oly. GOne gad. feo ite sbsteds «2 6 os 3s 3 
ers sak Directed Line Method eooeoeeeeeeeeeeeeeeee ee ee ee ee ‘i, 
bie POG ON ch wea eau Meth OG merciless sc ciettre cctehe aa hehe eiseenet ciate S, 
Aha eee} Sum Of Angles Method eoeee@e#e3#e+oeeeeeeeeeeeeeeeeeeee nee 
IAL aA Line Crossing Method @eoeeeese#oeseeeneeseeneseeeeVeveeeeeeeeeee 5 
IEE TN Pee Olyoonmand ee itswowe thede tmnt ton me cee a aan | ees 
did i eee Re Chiatn Mono Cone swith tie y-axis foee ce see ec 29 


1DoN os Sia) NOMS EONY OXmEO ly SOllmrsuster cre otetsccte citrere ere e cisterns ar ere ete mo 
ludk ibe eie) Pie@assOCtarccOmMinimatmregularmGraplim se cee ee ee temno 2 
Brie sec The Associated Complete Set of Monotone Chain .. 36 
III.4 Binaby Seacch #irec ot Monotone: Chains) .0..65 <6 or 
PELs) iMepleine|WMusheaiejn) he (Giaciels) 164) eletel (Esra ase oonaess 4 ay 
lit. oa Polygon enclosed=in a Triangular Region’... ..66. 44 


Pit. 6b Regularization of a Polygon relative 


Omame Lc NO Cuber cts tone falc) sists lot aisiereretctels (oolc ter cist ctelerenet ens ~ 44 
Pie Cc Miata git tig SUDO Villsd OlnwlmmiMetsteletete ists oisteleletsters-ctciciciceeaen eS 
iio NerohbouLhOogSmOrny CRt Cosme etctcter.tsts tits lel sis He dsnooa we 
LiL eoec Triangular Subdivision CL) Siatrre techs tt NN ge oct ce ead G 
JEW IE Slee Triangular Subdivision CAE Gane A anh h ae eae ~ 47 
Pit. og Triangular Subdivision AGERE werent yer A Pe gic ate cunts tees, 


Rite ol Tebvanguiac subdivision (Tse) aap teemerelst teres coher S 
PL. 6 1. Final Triangular Subdivision Tu 9) Ae Mie. opie eee CnC 
Pree, Hverarchicalsiree Structure of Triangles) ..-....6 49 
HIP ENE teks) BasvemsectuonPReGtang)] Gse—slevVe las sinc. sieves sso seo 
iieleeo Sect lonerectang lese =a bev Gler2) ccreretsyeielw selelets pc lcicis OO 
Eel]: S@GCEMOME REC tang Lose — CV Clu sie siere satatstaldinie s) «ecole kes EO 
1 AIE are Séction Rectang les=—) LEVel=a4 nce. sc cece seis ess 0 
LII.se Final Enclosing Section Rectangle —- Level 5 .... 5/7 


Teil HierarchtcalesStLuctures Of sections Rectangic® <<... sow 


Wy Le 


> 


. 


weer ese see es Steer aeas 


7 ' 
a hy ~™ _ rar © 


> acsoees SlQneEs? To Sie5029F 295T Iepitstessit | f.2EE 


=> _ 
PU ages < a, SS Be 


“4 ia Phy ye aie naa ae Bere hd = foe ba 


gen ules dbus pb eaidete.ghe ss 5S, 


oupmn see acsedeneenen etter rs “tet 
snbulvin tsa oes y) MASUR eb igno2 <4: | tute moeylee a 


boweeébantiacen ELSE : iste apnseaem are? & Pyper 7 
othe tas Be Aetae ea cae: “serpy ted Katnct—iee ttt 
wieneesvess $QE92 f6fupod Inakart Oeieidoens at? ‘OEs t 
.. afekhd adgiidet 76°98@ =scfuhtO bevekvoest eff DE. TEED 
see dead eee BR 2EBS.onGITIOM to asi? dowent yusenla “®.21t 
Loeece enbede ees Get Oe arian). =o octialoptsiat @.50t 
(iedcca. SOPREE Gettpoets? ant tsecidas nmpyta? e8.2ER 


seigg¢let nwivioy « to sofyssicelinweA dd. 3EE 
Seleee SS See eae pe 86 oe CE 6 eee es Fes OSS sitpae23T | os 


Set eeesebesderSenetey TD NOTeeusbews Isiegiela? oa. Ft 
capbcwa tise pedans~0s~ BOOED SEY 20 chetdagedtyie = 80, TE" 
064068908 F488 04 We Ae ri HOLGINSPAlR 46lupaetyt 23.25% a 
wtetavepaiesaccas- + Si watein Pie yetgensive  20,RE 

Tere tee ee eee ee ee (Ee) > eOis:vibdee® sealepaassT pa.rty : 
oécseakseasscesacas FP nGleieibies aelugaeier 86a EE 
wseheegeddaces, WD toilertndod tsiupasiet! ieadd . 26. 35E) 


ive vtetoawive, Cf loved - aaipisa ven nolgnse. 2:88 08. IIT 
O46 600 Deedes a> se Bde ¢ éovet — 89 <7" 
a@ ree ateraeteete ae ee £ Leva. + 


ee ee & 


dere @ ~~ * 
stew 4 


Figure Page 


Piel. Ode ecOncartmment, of POInt Goin more than One secrion 
Rectangle eoeeeeeee#eeeeeeeees#eeeeeeees5s#s:keeee#e#eeee#e#e¢@e 59 


DU tODBmolitnecmOl theshierarchical Structure wc 6 sslsccs Do 
jE Sh PeWOLs tecaces POLYGON mene mC.) GW )e ictus icusucicleteleistsielelarsm 0 = 
~eetnGo ae Tas 
ee? Dee OLY GOnat mA DO LOM ING ELOnmae cteleterstarcrecetsts clots ciel ste sie: tie enmen 07 


ie a A Given Polygon eee7e5s3e#ueeeecgretweenreeeneeesee#ee#eee 


Ijivil2ze sAnnulus Approximation of; Polygon’ P 
C250 BMECO Set OCEAN CE ic eteieate sc sreletelsislelc ese state etatere co 


Tit. 2a) Annulus Approxamarion, Of Polygon P 
(SO0SMELeErS BCOVEranCe.) mci oie orci tis a0 e scscel suai s, bas he ei oiel sie 


WIL GILSA A Set of Sip los eeeseeeeeeee#e#ee#ee#5egegee8#eeeseeeeeeeeeee e 76 
ie mn O Cia tao lmao m NOUCMO reram > Cl) Dietciets iets ice sieletstorerete sis sumed) 
Tibyt socambatasSUruceurestOnbad SLEUDSTreew es sc: . ote sls cliche ceisler 


RAL Tb 15 StLip Tree Constenuceon eeees+eeeseeeeeseeseeeeeeeeese@ 80 


why 


i . - | 7 


7 
FO 


c 


= 


: 


- a a | _ P 
‘eo ere 
_ oe pe 
“Le - 
pokises sna estos ht nde lag sileiecans 
Reh add dale seaneth sh Odekuk ss shabgerees elpusd008 
id ate@e@ea@rver8 axusouxss feoidotetasii si° 70 gg32 au2 GO. a 
aes (61D $d) nopgiet #283 teaow A tr 
aa cin teonwhedemewtasnceemy eae PAREBOe agaviso A att. 172 
Jo dec > acl tnne sl y9aeh Seropylot agi.m7r 


’ . 
. i : 


ca ee0ee eee eoeaeeve 


‘2 e#6@466969°70¢e28¢00 RP Bes 


S noggrat 76 aoizariyongqéA sulsesA ost. ttt 
Lu] caekweanes vel Seu seOs Syren mm 2%) 7 
o»@e ee © ce 2 @# 61% _ 
{ aopyigd-t6 adi cen xo1ge A suiwonA SSI.3ER> 
.e'G@@heee@re osa2@4 04446 se 6 (3° hs 41a sinzen O€) - -_ 
a es @e@=«@ iewreeenwr ee & . 4.0 2 oe. + & m= &e 24 78 +c 73a | sti, rn 
i ‘ojmees eebecceecsces GE322 6 to Shon 6 to 3emx0t GEL. i 7 
er pcnnvae Sob fate eotet 8ttouIe asad ats Lit 
noitoysssnc) souk 4722 bt.312 


- 


> 
- 


CHAPTER Rar. 


INTRODUCTION 


1.1 Background 

The. concept,,of a geographic information. system _jis 
essentially that of storage and retrieval of information 
aDOuc  Specitic Jocations Orsregions,,uSually. on the vearth"s 
surface. In such a system, descriptive data (such as 
NouSsingminftoORmMalwon,srainrall, sete. juts structured according 
to these geographic locations or regions. Because of this 
structure, rapid retrievals of descriptive information about 
a geographic location or region are made possible. 

The boundary of a geographic region in geographic 
information systems is frequently represented discretely by 
a finite sequence of points. By drawing imaginary cords 
between successive points, the resultant polygon serves as 
an approximation to the geographic region. The discrete 
boundary points are frequently obtained from their source 
(e.g. Maps). by. means of computerized digitizers; and in 
PEac eco, tne mnumbeierOL sD DOiIntoesLCOULreC ma LOm mallemaCcunate 
approximation is large. 

Pnewstoring. Gescr IpllVe ss datamesaccoclat cds wi thee sole 
geographic point, it soften becomes necessary .to identity a 
Bec Ones (DOlyGon )mewhiCchwecOn tad Nomeuhat pO tr. ThegeCcOstas0b 
resolving the identification problem, more typically called 


the point-in-polygon problem, dominates the cost of data 
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sivameih aa> . (Weber. Sittgeveper) cil’ a7” nobremixorgge iis || 
eagvec 36404 furs? it Sebsittes yitadarne |! ov attflod . ytefimedd 
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entry in most geographic information systems. Tce see thus 
problemithateistaddrésséed in this tthesis* 

The point-in-polygon problem arises in a variety of 
other applications including such important ones as image 
processing and pattern recognition. Whereas the size of the 
polygons (1i.e., the number of sides or vertices) may be 
Small win ssome™ Of) theses applications, in this, thesvs) the 


Eocus esha bl bevionwlarge ipolygons h(erg 281000 isrdes)i: 


1.2 Statement of the Problem 

There are a large number of computer representations 
for Lreqrons@andmpolycqonspilid, 155 1 7 flee 222 4 2G 28 esos ols 
Grids (raster), curves, points and polygonal lines are some 
OGM the? forms #useds Fommthempunpose voterthirvse fithesiss;tiwe 
choose to define a polygon as follows and in so doing we 
Testractenthes=sSCODe™mOLmthis sthesise ton ethoses applications 


Wherter thls caetinition dWSuvalid. 


DEF'N: Rete P=((\Vany Vogts? nt 1) eOcmamcedUCnCe Of DOlnts min 
2- dimensional Euclidean space E 2, where Vn+1]=V1- Denote 
the coordinates of Vj by Vj=(xji,yi) and let (Vj, Vji+ 1) be 
thes sliness segment, —jJOlningmeehnem point. V4 sand SVi7iem for 
iat, pic ipiete veteny ile Then P describes a polygon if OV Seal) 


does not intersect (VeE, VeE4+]) whenever 1 + k. 
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DE aN m oe DOLLY GOns DV jure) Vinod lS eCOnVex alt Ors any. 
two points contained in P, the line segment joining these 


two points is contained in P. 


DERN -sepAeceSty pOintyZiis fasset, of. coordinates (x;y) “in the 


cartesian coordinate system. 


DEEN wee tne point] 1n-polygong prob lLemmrssdelined gases ro! lows: 
Given the closed polygon P, determine whether or not the 


GeStmDOlnNtaZzeUs scOntalnedearh yr. 


For solving this problem, there are four possible 
relations between the given point and the region. The point 
can lie outside, inside, on the boundary, or coincide with a 


boundary point, of athesredvonvas gd | lustrated sine Fr ugune sisi 


oz 
(i) 


FIGURE I.l: Positions of a given Point 2=(x,y) 


Relative to a Polygon 
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AseiVlustrated in gFigure [.1, we say 
Case ls Joe Oo Ze is snOtscCOntalnea im regione.) 
Case il. ZeCeP (2 1S contained in’ region P) 
Cie cwi lence = Ze lomConralned ii ered lon ar) 


Case iv. 47 Ger (21S contained in’ region &) 


Ones Ccases that 2s 1s a boundary, point are distinguisned 
above because some algorithms described within give such 


cases special consideration. 


Poe Outline Of Thesis 

The purpose of this thesis is to survey methods for the 
solution of the point-in-polygon problem. However, no 
attempt is made to include in the discussion all works 
concerned with this problem. Instead, this thesis describes 
only those methods which are considered to be significantly 
different in their approach or performance. Since such 
considerations are necessarily biased, we apologize to those 
authors whose work we may not have referenced. 

In the description of all methods herein the intention 
1s always to present, only those “notions “essential to 
understanding the essence of the material. For details and 
for the handling of exceptional cases, the reader may need 


GComne SeretOethesOLlgindlesouLce. 
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Chapter II reviews direct methods where no initial 
PECCeESSIngmeOL ss thes ssidess= of, sthnes polygon | iss erequired, 
Extensions to some of the algorithms are discussed. 
Typically, these methods require O(n) operations to resolve 
the point-in-polygon problem. 

Chapter III deals with methods that require pre- 
processing of the polygon in order to speed up the search 
times) Typically fthese)methodstirequire 10(n);0(log2n), sor 
OGlogan n) operations ¥rtosferesolventethe st * point—-in—polygon 
problem,but in exchange they require the construction and 
the storages ofwasy more -complex® data istructures 

Chapter IV reviews methods which quickly eliminate the 
areas which do not contain the given point. Examples are 
presented and implementation factors are considered. 

Chapter V summarizes the previous chapters. 
Conclusions’ tare midvawns and erecommendatronssytor pofunther 


research proposed. 
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CHAPTER II 


DIRECT METHODS 


2s ee i ncCroducvion 

This chapter describes those algorithms, for resolving 
the point-in-polygon problem, which can be applied directly 
to the polygon represented as a sequence of the oriented 
VEDEICES@EV] Fics ty Vy, eescer Chapter a igr THeSeemeLhods sare —to 
be distinguished from those of Chapter III which require the 
polygon first to be transformed into different represen- 
tations before the methods are applied. 

The methods described in this chapter are therefore 
applicable to™athose applications where a polygon is 
subjected to frequent changes. No preprocessing of the 
polygon representation is required before the 
point-in-polygon algorithms of this chapter are used. 

Two of the methods are restricted to convex polygons. 
Polygonal regions which are non-convex, however, can be 
SOC mincOmeDe .Un1On OL sconvex* polygons ane! Gr. CO) NOrdbeck 
and Rystedt[23] where a non-convex polygon is made convex by 
addingetOilt an appropriate number of jtraangles., “Ing those 
cases where many new triangles must be added, the two 


methods affected may become very slow. 
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2.2 The Directed Line Method 

This method, described by Aldred[1] and by Nordbeck et 
al[23], may be applied to the point-in-polygon problem only 
if the polygon is convex. 

Proceeding sequentially around the polygon in a 
counter- clockwise direction, at the i'th step, the method 
determines whether or not the given point Z=(x,y) lies to 
tremletti¢gor tomthesright ofithes line, directed, throughtgthe 
ventroes ,Vj=(xXj,yi) and Vi4j1=(%j41,Yi+]1) of the polygon P 


(see Figure II.1l). 


Vig Gap Via) 


BLIGURE LLL: Directed Line Method 
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comparison). 
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Signiticantly less, "since ‘the’ ’method* "terminates with the 
PLGStmOCCURLCNCe: Oba i (6X7 Vi) esa. 
2.3 The Area of the Polygon Method 
The area of the polygon method described by Aldred[1], 
like the directed; line method, }is restricted to polygons 


that are convex. 
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BPhoceeding <sequentialilly around®’*the polygon! in’ a 
counters tclockwiselidirection; atthe Sith step, athe method 
determines the area Aj; of the triangle determined by the 
given point Z=(x,y) and the two vertices V4=(xj,y;,) and 


Vi+i=(Xi¢1,7Vit+?) Of ythe wpolygqoner «(see4Rigures«DI $2). 
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FIGURE II.2: Polygonal Area Method 


The area Aj is expressed using the determinant notation as 
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The determinant is to be interpreted as pairwise 
multiplication of elements in the two rows and_ then 
subsequent additions as follows: 


2Ai = [(xyitxiyititxit1ly) - (yxityixitityi+1*)! 


[x(Yi-yit1)+xi (yiti-y)+xis1(y-yi)] 


It is easy to show that the area A of the polygon P 


satisfies 


n 
he eo 
i=l 


The summation above may be simplified as follows (but note 
that. thiss simplifrcarlongpis @nol ‘germane to. the point—-in- 


polygon problem): 
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cigs Ae On BSpepgococ pity We Va 1S IF ee 4) O(P), and since 
all of the Aj's are not of the same sign if Z ¢ P, it follows 
that 
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Thescosteom Aldred's .methoa 16 


Cy =Fn(3sCm + /Caar Ca) 


al 


(where Cris “the “costwoft multiplication, C3 is the cost of 


addition or subtraction, and Coe is cost of comparison). 


As an extension to the above method, it is clear that 


2 CBP #1 Ee AME er l= 1, 24. tee al 
Ze CyO Pit eAye = 0; cOmesome i), 


ZU SPSL AS <0 ror somes: 


Hence, the cost of the method is reduced to 
C= n(3Cm at ler) an Ce) 


inethnes worst case. 


2.4 The Sum of Angles Method 


This method, due to Nordbeck and Rystedts[23] 


proceeds 


sequentially around the polygon, determining at the i'th step 


the angle 9; subtended by lines from the points Vi = Can) 


anda te — (X41, Vie) sto: the point, Z=(x,y). (See Figure 
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BILGURE I1.3:> )sum of Angless Method 


Toewangle 07 91s given by 
6; = tan-l(yisi-y)/(xi¢i-x) - tan7l(yq-y)/(xi-x), 
where -1809 < 6; < 1809. 
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determined that Z ¢ O(P), however, the effect of round-off 
errors is not otherwise critical to the ago thmeethats 1s: 
Z€P if @ is "close" to 3609, and Z ¢ P if @ is "close" to 


00, 


A considerable improvement to the methods has_' been 
Suggested by Hall[{14]. Rather than accumulating the partial 


k k 
sum S, 07, the Sine and cosine: of ye 6; ares saccumulated 
i=l i=l 


msteGadquethen, at Ehervertex Vien, 
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The improvement in this approach arises from the 
Observation  thatercin (0; a)emand)=COs( Cr+] )memMayan pe =. compuced 


using vector dot and cross products. 
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Ene eVeCeLOtmcdo le prodlctasat is ties 


(VErZ). (VietanZ) = |Vk-2 


Vk4+172Z | .Cos O41]. 


Piepvector cross joroducthesaticties 


(V_-Z)X(V_41-Z) = |Vx-2 


Vig 4—Z |.Sin O41]. 
Ondysateethe, finals step does 4 trigonometric \routine, need: to 


be invoked in order to compute 
in 
@ = tan-l(Sin 6/Cos 6) , where @= Dy 
i=l 


The «computation of sin( 64); Fcos{ } 67 Jdiv.en Sin ( Yen, 
1= 1= l= 


is 
SEEN SS On) requires Poe enero licat TONS, a Umm aGdite tons me 
1=] 


Givisions and lgsquare root evaluation. § Ignoring the cost of 
thew fina Me@sean--» computation, the cosy, of Haldi'ss method |ic 
given by 

Gi=en( 13 Cot elOCa ts 2C Get aCe.) 
(wheres Cy 1s “tne cost ‘of multiplication, Gs, 1S (the cost. of 
SBddi ton Or ssubEraction, Cami Sethom cOoStseole division, sands Ca 


is the cost of a Square root evaluation). 


Deol nesGrossing Method 

The method has been described from a variety of 
Derspectives bysa mumber Oolsdi Eferents authors) [1273747570777 
Oslo Oy 22 po 2 lus It can be applied directly to polygonal 
regions which are convex, non-convex, and even to those which 


are not simply connected. 
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Basically; mthe, methodsconsistssimirstilol, passing)ta, lane 
through jthesyspoint? Z=(xpy)aeandaaparallel” toPfithes@rx—axis+ 
(Al thoughehanvyiaelhincarthrougheeme wiliae doyivasnhormzontal fier 
vert icalggline tis) tcomputationalily ebest.. The number of 
intersection points of this line with the polygon boundaries 


is then tabulated (see Figure II.4). 


FIGURE II.4: Line-Crossing Method 
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ConputatonaDlyyiiat thee ivsthiisitepal 1=12 ae. ae the 
sign of (yi+1-y) is determined and compared with sign (yj - 
Ve)re If the signs are different, the polygon boundary has 
likely crossed the given line, and x-coordinate of Z is then 
Compares wltheetne) x-COOrdinates “Of §Vi,. Via) Lm sOrders to 
determines thes pos turon @ofivaerelative to §the “seqment, Vj, 
Vier] #e loa Ze ise Contained between thew extreme sx-coordinates 
Ol V i pmevn tie) then thes vertices V7 and V7.7 are -linearly 
interpolated mitosmdetermincemrt SiZeeaeseron@icr lett of the 
polygonal segment’ 9(Vj,Vj4+1)- Special attention must be 


Givenmcomthe case that the sign of (yi4n—-y) eis 0. 


The “aligorpthmeiny greater «detail 1s ‘given below. The 
algorithm returns the value CONTAIN. If CONTAIN = 1, then Z@ 
is contained /in the polygon; if CONTAIN = 0, then Z lies on 
the, boundanyeokethe polygon; Wand finally Atif VCONTAIN = =], 


then Z lies outside the polygon. 
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een 
(Comment: 
go to step 4 
else 
(Comment: 
go" toustepngo- 
he SCARCE lh OA De) 
then 
(Comment: 
set CONTAIN = -CONTAIN 
go to step 9 
SOmtCO jst epee 
bimax(%7 17 x) 
then 
(Comment: 
gomtomctepsg 


else go to step 6 


Zmlves left of segment 


Z diesiright tofmsegment 


17. 


>] Sandmsputeecordmthe 


top the vertices Vjo7 


Vertical cross-over found) 


Process next point) 


(Viena) 


(Vi-1,/Vi) 


- ee 2 
= 7 bet = i 
feaat 38 ele = oe ‘ 

| ee 

er a 

(pkuod. tavo~aebils Lanisye) =Sasanso) capt 7 


> qase et Oe © 4 
eagle .« 
7 
(ieion dean 2363025 -3ne~aes) 7 ' 
| ~@) G84e 2) Pp 
(cx) cast) aie 5 e ah 
eee Be) 
© 
HpPepaghl; titeess Bo sted atts taqaeneas = 
WEATID- = OEZTWIO tec a 


& Gose"G? OD 
2 gate GF op wake 
ae ee 
nde 
(g¥ad-2¥) Goumpse 9p siya 20ld 2 sdieames o 
| eget ste 5 6 
a qade atop epis =e 


7 


18 


on (Comment: Z lies in rectangle containing segment 
(Vi-1,Vi), and therefore interpolate) 
eH Sl Gur Ore ly s(xie— xy iy) eyes CXR) Pty x— xe 
Tf (S> = —1)"then“set®F =*=F 
thee Lev rs=—'0 
then 
(Comment := Ze 1lies= ony segment= (Vj2i,"9V4) 
return (CONTAIN = 0) 
elsewgorto step ss 
Bi Loser) 
then 
(Comment: Z2 lies left of segment (Vj-1, Vi) 
set CONTAIN = -CONTAIN 
on ia ob Olan ee) 
then go/ito step 2 


else return (CONTAIN). 


Since interpolation is required infrequently (at most 
twice, sib. P 9 is ~convex), “the “costs of) the) method = is 
approximately 

Ge" Cae Ca)) 
(where Cy, is the cost of an addition or Subtraction, Co 1s 


the cost of a comparison). 
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2.6 Summary 

In this chapter, we have discussed four methods of 
determining whether or not a point lies in a polygon. These 
Mmethodsm, arewasumarrzed ein Biterms “Otss theireyworstwemcase 


computation cost and» in terms of restrictions in Table I. 


Method Point Processing Cost Restrictions 


Directed Line nN (3Cm+5CatCc) Convex Polygons 


Area Of Polygon n(3Cmt+5CatCe) Convex Polygons 
Sum of Angles n(13Cpt1l0Cat+2CgtCoc) None 


Line Crossing n(Ca73Co) None 


TABLE 1: Summary of Direct Methods 


It is clear from the table that the Line Crossing method, 
in the worst case, when the point is within the polygon, is 
ERee Dest eon all. the we Direct, Methods". discussed, ied cole pbenverel 
itis not restricted to convex polygons as are the Lirst «two 
methods and 2t. 1s Significantly faster than the third method: 

When a point is outside the polygon, however, each of the 
first three methods may terminate early, whereas the line 
Crossing mimetnoOds s COntInucsm tOss DrOCcesss fall aepOlncSes Oleic 
polygon. Thus, 1£ the majority of points™ tested happen to lie 
outside the polygon, it may appear that any of the first three 


methods may be superior, on the average, to the line crossing 
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method. However, aS we shall see in Chapter IV, initial 
inexpensive tests can be applied to resolve the point-in- 
DOlycOnmeproblems, Lor “most wpoints which lie outside “a 


polygon. The superiority of the line crossing method, 


therefore will prevail in all cases. 
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CHAPTER@LDE 


PRE-CONDITIONED METHODS 


Spe A eLOcuction 

All the methods described in the previous’ chapter 
require O(n) operations to resolve the point-in-polygon 
problem. Asymptotically, better methods are available, and 
these methods are the subject of this chapter. 

DRevLemwis  sOwevVer ma pLice tOmDaymiOnmthicmopecd mal 
the methods require that the polygon be represented by a 


Certain, Often’ complex structure. The storage requirements 


POr this Structure are -orten’ exorbitant (O(n2)) ANGTCOs ESmoOr 
achiev nGeeathiis, structure (preprocessing) can be high 
(O(n2log 161g) i bee Therefore, these methods can be practical 


only in a static environment, that is, in an environment 
where the polygonal boundaries seldom change. 

There is one additional drawback with these 
asymoOcoucicaldy = siast we metnods. Because” of" "the “complex 
St LUCLULCmDy = will ci d= "DOlyYGOns [US represented, = accCessm and 
Other soverhnead costs reduce =the "practicabivity “ort "fast 
Metnods= LO problems "whereml is**Vvarge =) “Inem threshold tor 
n, aS n increases, wherein fast methods become superior to 
classical methods (e.g., the line-crossing method) is an 
important implementation consideration. This matter, 
however, is beyond the scope of the current investigation, 
and is not addressed in the following pages. 
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Five methods, namely, Swath, Chain, Triangle, 
Rectangle, and Strip, are described in the following 
sections. 

Bach, Of theysectiong; includes ‘a “description ‘of the 
required PEEPpECcessing Of) <a “polygon, followed by a 
description of the point-in-polygon algorithm. Costs and 


storage requirements are also given. 


3.2 SWATH Method 

aes method is due to Salomon[29,30] and ee 
Shamos[3l1]. The preprocessing step proceeds as follows: 
Preprocessing 

SWathS ©:dLeG= tOGbMeG er Dys —CONSCRUCUING @horvzontal lines 
throughy each vertex of. the polygon, The swaths are then 
ordered by-.decreasing y-coordinates. 

To illustrate the formation of swaths, consider the 
eight-sided polygon given in Figure III.1. This polygon is 
described by six swaths. SWATH(1) contains edges (5) and 
GGpeewithinenchemy—InCeLValele lo. 7ue2. 2) peONA LH U2 econ tains 
Sages oS ye 4 )meands (6) gwithinethe y—intervale (2.2 7.9) ead 
Sinai ye SWATH (6) a .cOntatnsmscages (4) anda.) within the 
y-interval [-2.6,-5]2. Having obtained the swath 
coordinates, the swaths are ordered. This ordering is 


possible since the swaths are mutually exclusive. 


1 yy swath is defined as the region between two such 
horizontal lines. 
- SWATH(i) denotes the i-th swath. 
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Finally, the polygon segments within each swath are ordered 
from left to right and indexed. 
The complete swath method is summarized by the 
following two algorithms. 
Preprocessing Algorithm 
i. SOnte the@iwertices Mair (xirYi)s LS dy soste se sree 
according to the decreasing values of y-coordinates 
VolemnDenovemthetorndercdwlist. by ¥.(1)), Ml=1eia ace me 
from which all equal entries have been deleted so 
Gita Gameoan. 
aa) DEnoLew@es cache consecutive! By tapalr sioL . endpoints 
(vcGis) py Gite 1) Ame SOL eel =l) 25. rey Mody (with) =Y¥(1) 
Obit) ae yar WAC EE (ela) rs 
game Within SWATH (i), establish a Teéfit-—to-right index of 
all line segments which intersect it and store the 
Pidex in SWATH ait ly)epae jy 2 peek ets 
(comment: SWATH(i,l) is the number of segments in 


the SWATH(i)). 


Point-in-polygon Algorithm 

The algorithm accepts the matrix SWATH(i,j) and a test 
DGIncez=\ xy) Sasmanput;s andM@meturnsesas foutpuGeethe= value 
Lirire! aeit eZ is) Scontained = indethetepol ygony erandeevea lise. 
otherwise. 


3 SWATH(i,j+1) denotes the jth line segment within the ith 
swath. 
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imeecOMpare y against.yY(1)+ i=l). .0.,metosdetermine the 
SWATH (i)eecOntainanga ly Lf gno gswathyecontaimsie ys 
then return the value 'false'; otherwise go to step 
oe 

2-{),SCan, line; segments jinmiSWATH(i)esaccording to the 
COrderey dictated webby sesSWwATH (ae) ) jee 3 yest Let 
SWATH(i,Index) reference the first line segment 
wha chuyl i espeonwethes righth OfanZ. iy no @Ssuchweindex 
exists, then return the value 'false'; otherwise go 
LOOMS tepes. 

Steed felndex  —) smodulo. 2 ,mthen=return. Le ouew. 


ELsesreturn ~# false, 


Analysis 
[Preprocessing] 

ASSuMings Ny siethe numbers Ofa Sides sol the “original 
DOlYVdGH, wthecewwares ate most gnaleeswaths. Thusea(n glogmn) 
operations4 are sufficient to sort the y-coordinates in the 
decreasing order for defining the swaths. Within each swath 
there are at most n segments..(2.if the polygon is convex). 
Th theresbore. takesaat emost (ni.log@n) moperationss fomeche 
left-to-right ordering of segments within each swath. Since 
there can be O(n) Swaths, the total cost of preprocessing is 


O(n2log n) operations 


4 log n denotes the logarithm of n to the base 2 
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[Storage] 

From the above algorithm, there are at most n-1l swaths 
to be stored away. Int*each® Sswath7'@ there ™varertat@imost Gn 
segments. Subsequently, the total computer storage required 


by the (n by n) swath matrix is O(n2) units. 


[Search] 

Finding the appropriate swath containing Z can be done 
in O(log n) steps (using binary search techniques on n 
elements). In the selected swath, there are two 
possibilities. 

Prey the =polygonWis convex, fourtadditionaleicompar isons 
are sufficient to determine the position of the point with 
respects to ‘thes quadrangle™ formeda'by=* stwo. Ifine segments and 
the swath. The totalesearchVanetthrs #caseeusit( log pnestit) 
Operations and the cost is therefore O(log n). 

Pie Sthe= polygon 1Ss% notieconvex jpezenis: Compared) against 
each line segment in the swath until a segment on the right 
of Z is found. This is essentially the line-crossing method 
(described in section 2.5)* applied to’ the pollygonal portions 
in the swaths. Since these polygonal portions can have as 
Manys'as n-li sides’ in theeworst case,*the) cost) of this step 
ti =Salomon cee method Ganda Shamos! ¥methodimtcan” attarninva 
complexity MOfPeO(n)™ Operations! (however, G@sec refinements 


below). Bulerin a practice mathew numMbper Or sl INneCumscoments 


Sevéupsy sperode serngadD: [ado 
,esene (8n)0 2: ei dean tvewe (a y@ a) ado yd 


> 


souh ad nea 5 painiesiwon dime ateragoesgge onJ er Lori 
A ac seupindtse? dogese Vidal’ gated) eqste Ue goljo al 
aes  es6 S194) . ewe ‘‘bedssiee. wit int .(z2nenels 
.20i2i Citar — 
Smoeizagmon fedtohlitte WwHs «xe? =f newytoa aay if. — 
ftie gntoy ses To aeris cog sid onits7eb oF Insiaiiive a 
hms ednonpse sail “Gwi Vd becriga’ slodeybdup sit. cr saegane 
() + fh woh) ef Sése0 git? mi dowase [e203 ST . sige odd 
th pobji svote1sht 41 J2eo sar ‘bie E00 £96 Tao 
seninks PosnuwoD, at. T oy Renn” jca, =! aoey fog see 27° oe 
gipia. 2% oo myompee 6 Iifau deawe sala nl-tepense anbt aaae 
hodjett GateeSss—4ealT oef vilectrsces as 2iay Fito? ern Ro o 
efibbisoy fehepyiong. ad2c» beilqaa:(i.f oorjosaqatt aiden 
Sh Sia #eo enoidyor. lenopglog sass s9028  .shdewE aia & 
gintt Jo S205'dde.,, 2809 faiow Sit nk dsbie “<n es youn | 
iette neo bortsm 'gamalla bile beitom @'nomod 
ednemeniies 292 Zeveued) anoitarsgg. “tnyo) ae wy 
: 
meee rr, weds we cd 


atts Pn 


2 


in a swath will be much smaller than n-1 (2.0. two-site the 
polygon is convex), and the cost of the Salomon[{29,30] and 
Shamos[31] algorithms will be dominated by the cost of 
searching for the appropriate swath containing Lie 
Summarizing the cost of Salomon and Shamos algorithms is 
O(n) Operations in the worst case, but O(log n) operations 


for most polygons. 


Comments - Refinement 

An obvious refinement of Salomon and Shamos algorithm 
is Vavailable* Having determined the swath containing Z in 
O(log n) operations uSing binary search techniques, the same 
techniques can be applied to the ordered line segments 
within the swath to determine that first segment lying on 
the right side of Z@. Since both steps are now of complexity 
O(log n)sisthe @cefinedmalgorithmrtvs: iof S&the «same. sconplexity. 
With this refinement, the algorithm becomes essentially that 
Ge DODKIn sand «Lipton ss] % 

Nota, Showever,:sthat™ Dobkim Sand (Mipton | presentestherr, 
method in a more general setting. The problem considered by 
cChemeiceethaty Of identilying ss that = polygons (11 s.any)eewolch 
contains: a given point from¥the setmofr- all “polygons *formed 
pyetthes intersections: of nn Slinesi int apeplane: But any 


prescribed polygon is determined by a set of ordered lines 
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and the intersections of each line with its predecessor and 
Successors only. Thus, the point-in-polygon problem is a 


special case of the problem considered by Dobkin and Lipton. 


3.3 CHAIN Method 

This method, discussed in terms of planar graphs, is due 
toy Lee and Preparatal[is]| who # utilize “certain “geometric 
properties of a chain for solving the point-in-polygon 
problem. The concepts underlying this method lend themselves 
to many different applications. For example, their method can 
be used to determine if a given point is contained in a convex 
region in 3-dimensional space (see [18,20,27]). 

Before describing this method, some preliminary 


definitions are required. 


DEF 'N: Given a polygon P with vertices Vj, 1=1,2,.....,n, a 
Chain Bwitierespect tomP)etis detinedstol beiman ordered isupset 
(MO iene = = -)7Vi 8 Offa Vertices 408 0P- InplTClit winch s 
definition is the existence of an edge Cin a Joining the 
vertices Vi and ion 

DEF'N: A chain is monotone with respect to the y-axis (or an 
aEDLeralysee ine wm )meel cit nem Orthogonal DLOGeCelONS moO tacic 


vertices of the chain onto the y-axis (i.e. the y coordinate) 


are ordered. 
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An example of a monotone chain is illustrated in Figure 


Tyds Tet 2g 


RIGUREsi tied: yAyChain. Monotoneswith the y-axis 


DEEN ee cuee se tarC—\(C], Co). Ck Ole ChaimomwithierespectasLo 


the polygonwgP) is sicomplete, wishpte,satisfiesp,ithe stollowing 

COndit1ons: 

i Each edge of the polygon P belongs to at least one 
Gha lnsCy 

Ze Bachmchain Gy, l=, aye) Ky lon MONOLONGCw WILL respects LO 
the y-axis; 

3. Theschains) can be ordered such) that 6]<C 90.0. Ci where 
Ci<Cj means Cj lies on the left of Cj. 


Thise caefinition ~ensures’ that. the edges or Cy > do notwcross 


those of Cy, i+), and that. a total ordering) of @thegchains is 


Maintained. 
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Lee and Preparata's method begins by constructing for a 
polygon a complete set of monotone chains (the preprocessing 
step). This construction then permits «the fast’ resolution 
of the point-in-polygon problem. 

Only convex polygons lend themselves naturally to the 
construction of a complete set of monotone chains. For 
Other polygons, 1t is convenientw Eirst \to transform the 


polygonPinto a Suitable reguilatmetorm. 


DEF'N: Let (Visi 6.eV,-) Dew de SOrting: OL@ethes vertices 
(Vile aw ny) Ob ma polygon Piso thatyy (Vi Sy (Vo. Sy (Vico 
A vertex tia is said to be regular if there are integers 
i<j<k such that edges (Vai ve) and (V5*, Vie) aren ab sone & 
ell of ethe vertices with «the exception jor Vi~ and Vy;,> are 


regular then #2 1s sardeto be, a regularsgqraph. 


For example, the vertex H of the illustrated polygon in 
FIgQube lili, oa slssnot cregqulan «whereas Geis. 

The transformation of a polygon into a regular graph 
consists of drawing a minimum number of additional edges 
betweene vertices. Of 6h) sOwthnat, for cach Originals cdgemor oF 
there exists a monotone chain in the graph passing through 
V,* and Vp* and containing this edge. An example of a 
polygon and its associated regular graph are given in 


Eiroucess ivi, cavanid Llleol. 
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Non-Convex Polygon 
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FIGURE III.3b The Associated Minimal Regular Graph 


©). 


The regularization scheme’ is déscribed as follows: 

te Ore the vertices in non-decreasing order of 
VWeCOOLd ina Ges as Vii mi ome Votstety Vine 
(Comment: Vj* is the vertex being processed) 


2s Set j = 2 


Ss Set 1 =j-1 
key +> 1 
4. (Comment: Descending Pass) 


4.1 tt theremrsmaneecdgou(V,* 7) V7 4) 
then™=qo7to step 5 
etse see —me1e— 

2 a0 
then go to step 4.1 


else 


go to stepr 4% 3 
4.3 (Comment: Draw a line segment) 
Ie TAs v5") dOeCS NOU InNteEsectmanys OLlermscedment 
then 
(Comment: such a segment must exist for some i) 
draw segment (Vj*, Va) 
go to step 5 
else 
Seti = le 


go to step 4.3 
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(Comment: Ascending Pass) 


Syl 


If there is edge (V5*, V,E*) 
ehensgoyrto: 6 
else set k = k+ 1 
Dt ek <aen 
Gh enyiGOmto ysitep, Sys) 
else 
setake=ijntit 
Gowtowstep? 53 
(Comment: Draw a line segment) 
Ge (V5*, Ve*) does not intersect with any other 
segment 
then 
(Comment: such a segment exists for some i) 
draw segment (Vi*, V3*) 
GOetosstep 6 
else 
setekpes ky +h 


GO, tO FStepiaor 3 


See Gh Sage ar ab 


ike ye eae 


ENenedOmLOeS cel 


else (Comment: all vertices are regular) 


Stop. 


n> a4 32 S68 
L.@ qazge oF Op neve 


sels 


(zaomppeventiis Wasd +faamitgg! E.é 


(* @. 7) 21 


Yargo yre djiw InseTrssnt, 100 Be0b > La 
soonrse 7 
nist? 
ft amésa “62 afeixg tadtmse a iate <Fneeme) 
(Fi. .%,V) Shetipas wast 
a 9sste Ga oe 
sele 
oo eo _ 


iat, kg SOE Gd ae 


35 


Once the minimal regular graph for a polygon has been 
Constructed; all) uniquemmonotone = chains (Cj, C 97... Cl) 
passing through V;* and Vy;* are ordered so that Cj lies on 
Orson the left or Cc) Pitre iy <Q. Thesewchains form a scomplete 
set Of monotone chains. The monotone chains for the 
exemplified polygon are given in Figure III.3c. 

The complete set of monotone chains C can _. be 
represented hierarchically in a binary search tree. The 


algorithm for the binary tree is described as follows: 


Givens am cCOMplewe Set, Ofeechalns = Ce— m\ Cis, Co, mee Cm) Gad 
binary tree of m chains is constructed by selecting the kth 
element of the set ,(where k = |m/2H1 ) as the root of the 
frees These k= element partituons: sthe seu © Into sa slefte 
SUDELCe (Cl Meo. + -.ch 1) wen se lementomvess ethan Cry ands a 
GicghtesuptLee |(Ce lily. sy cm)  Wltae eclements= greater —tirang Ci. 
Pach subtree 1S )lurtherampartitioned by its, middie selement 
into smaller subtrees. This process is repeated until each 
subtree is null. For the examplegat hand,. thesearch =tree 
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in Pecheeagtree, the chains ~~ associated. with the left 
subtree of any node (father chain) all contain smaller 
indicess ande@thereforen liesonpthet left for) tha stathemichaint 
A similar observation can be made about any right subtree. 


The pleaves ’nofOel thes treet denotew@icontainment Fini {4)70 or 


exclusion from (—), the polygon. 


m4 
Cc, Cc. 
mal Be Cc. = 
re nue NE 


Chain: C, AHCFE 
C, AHIJGFE 
Cc, AKJGFE 
f, ABDE 
ABDCE 
Cc. ABCE 


FIGURE III1.4: Binary Search Tree of Monotone Chains 


38 


[Search procedure] 

The search algorithm begins by checking if the y 
coordinate Of saigiven) test spoint Z=(x;y)slies within the 
extreme y-coordinate range (y]*, yn*) of the given plane. 
TheePEOCceSss ingwOfethe pointazaismterminated sit 11s stoundmto 
be outside the defined region. Otherwise testing commences 
with the root node of the search tree and proceeds downward 
toward a leaf of the tree, whereby the point-in-polygon 
problem becomes resolved. 

Testing of a node at each step consists of determining 
Piece  pOlnte| Tesuonmthe lettwor sons thessuiGntmoOmethes chain 
associated with the node. It is interesting to observe that 
Ones Guminat von of the search algorithm,.i4 sun gue. palLmsor 
GConsecubives chamns icq and G7 sq. containing the. (est, polntaare 
determined. Note ethatyui fi ether pointe 7a se Strict lyeon. the 
left or right of the complete set of chains, then only one 
chain is identified. 

Discrimination of a point Z against a monotone chain 
Gur cmmOeLeLMini nga 1. simllesmon othe Slet ior ona themrignt sot 
pBhemchain)« requires Giust thes identiuticationgore the edges s0f 
the chain nearest to @. Since the chain is monotone with 
respect to y-axis, binary search techniques can be applied 
to the ordered y-coordinates of the vertices of the chain. 
Once the nearest edge to Z is identified, discrimination of 


Z against the edge is straight-forward. Tllustration of 
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this discrimination of Z against a monotone chain is given 
Hrirekroire LI 2s 

Analysis 

[Preprocessing] The P init val Mi processing™ or e*the eepolygon 
requires the sorting of the y-coordinates of the vertices, 
which takes O(n log n) operations. Scanning the edges of 
the" polygon” "in¥=*choosing™ monotones**chainsee takes |= O(n) 
additional operations. Lastly 7. the =-constructionevol=+the 
complete set of monotone chains requires O(n log n) more 
Operations. Hence, the= preprocessing™ cost, Aseat-=mosu Om 


fog=n)* (see Lee =and Prepdratal 18) tor -additional“detalis): 


[Storage] Since at most n chains are required to form a 
complete set and since each chain can have at most n edges, 
the total storage requirements are bounded by O(n2). Lee 
and Preparata, however, describe a method for the storage of 


selected edges only, which reduces storage costs to O(n). 


[Search] Suppose a complete set of chains C is given by 
GCG, Gory srt Cle) USinOMeetne se brnaryesscauChmetrecom= Lorene 
Conn letem'sete=oL @chains, sat, MOstew loge k= chauns==needs@ tbe 
processed to determine a unique pair of consecutive chains 
Goeand=Cyy pe (with 18k) ewhich*contatn™ the™ givens pornc ez. 
Het e tne projcctions Ors the Vertices of Saeciains ontomtie 


y-axis form the sequence Y(V1) -Y(V2),-+¢+rY(Vm)- We can 
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apply binary. search .techniques. on. the, ordered set of 
vertices in order to locate that edge of the chain Cj that 
borders’ Z. This stakes sO ClOgem) =timeeeA fixed numbermor 
arithmetic operations and a single comparison are sufficient 
EOe deLermine@ On ewhich  sside eOtestnewm edger the  polncmt tes: 
Thus; tne wentires search process, requires O0(10g kk.) log, m) 
operations. HOWGVeGre eel COlbe che —Carliers discussiloOnmeaLe | 1S 
clear thatw there lare, atsamosts n—l..chains vin eC. Tnmaddd cron: 
each chain in the complete set C has at most hn vertices. 
Therefore, at most O(log2n) operations are required by the 


search technique to solve the point-in-polygon problen. 


Note: For. convex polygons, the .complete. set o£ monotone 
chains consists merely of two chains. The cost of searching 


in this case reduces to O(log n) operations. 


pete clangle, Method 

Ag sekructente 0. (ne Ogun) sa.LoOnmhthm bor spantitloningesd 
DoOlyoOvmer into sd trang ata subd1 Vis TOneel Ss descr bedm by 
Garey et -adld 21. Theil cy aloOb1 thet Les ta COnSsELUCtS eons Ee 
polygon an associated regular graph, and then an associated 
ordered sequence of monotone chains (see section 3.3). Once 
this structure is available, the triangulation of the point 


set which determines any two consecutive chains proceeds as 


follows: 
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Tne SteLangulation algorlthnme1s given as) Lollows: 


New: 


Bete the “vertices oOf the polygon formed) = by two 


consecutive chains have k vertices: 


1s Sort the vertices into non-decreasing order. of 


YecOOrd inmates, andedenote stiemuby V1 w,55V9 "is een V io 


Ze Place=the two: vertices Vj* and Vo. 1ntoOva stack 
Set) = 2 
Ore (Comment: Denotes cer CULECKt CONNEC eo Oe tne stack by 


Wapu WO 9 ene mer 3) 
Seta) =f yer 1 


Setuv* = Ware 


4. LigVemlseadjacenie LO WT Oute NOt sLOeW 
Een 


(Comment: V* lies on the chain opposite from the 
chain associated with the stack) 

dcaw dlagonalsen Vit, WOM Voy Wo) irene clony. \yany qe) 

change the stack to Ws, v* 


go to step 3 
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Ss TE \V* is, adjacent) to’ Wycbut not to Wy 
then 
(Comment: V* lies on the same chain as the stack) 
Repeatwuntilwei=lmsor Sinternal angleatata wae isaeat 
least 180° 
addsdiagonala( V4 pWwa_ 1) 
delete*W; guromstack 
set i = i-l 
end of repeat. 
AdGMV ta COmtnemtop Ofstheus tack 
go tosgstep 3 
6. If v* is adjacent to W, and Wj 
then eaddeaditadonals @( Vise Wo; ete ete, CN ar Wa ae) 
stop. 
ThemappJication of the) striangudation-salgorithm on: two 
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Lipton and “Tarjan(19]}e use this LYlangulanesteructure sto 
descr bempetiic#ents algovithimsaston mcolvingmeamevariety sof 
classical planar problems, including the point-in-polygon 
problem. A simpler and more practical adaptation of their 
point-in-polygon algorithm is given by Kirkpatrick[16] and it 
TSecit Seworkethatem s Outlined winesthis seation. 

Kirkpatrick's algorithm requires that the polygon P be 
enclosed within some triangle T. The polygon P as well as 
ENeumsECtTONMeXterlO@r sto, Pe andmintenion elrels iret ineamscor a 
triangular subdivision using again the methods of Garey et 
elec An -example of a pollygonigand Yits” triangulang. sub— 
diVUStOneealees OLVene gin hLoOuULeSs «LLi.oae) tOmell le 6c. The 
triangles composing the subdivision of T are labelled 
nO eer rh) pai $2) een t 

Given, Tan yarbitrary ws point 92, the problem = now 1s “to 
determine which triangle p, (0) pel eeany eel eet hemcrlangu lat 
Subdivision or Tt econtains \ za. Ue YA IS As (0) ana if 7 {0) is 
PiLeuVOrmtOme,— tiene Zits eCONcCaIneCOmMI) sor OLIeCLWise 1 La tounOt. 

The problem of identifying that triangle p, (0), liwany, 
which contains Z is simplified by obtaining a new triangular 
subdivision 7 (1) = {p, (0) el Le, Pp Roma ic. ea hoe of T where 
nj <gno. This néwssubdivesion 7 (1) is obtained from 70) as 
follows: Let V be an internal vertex of 7 (0) , Thenas Vas 
incident with some number of triangles of t (0), the union of 
which is called the neighbourhood of V. The strategy is to 


select a maximum number of such mutually disjoint (except 
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FIGURE III.6a:a Polygon Enclosed in a Triangular Region 


PIGURE LIl.6b:.) Regularization Of Polygon Relatives toma 
Triangle 
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FIGURE III.6d: Neighborhoods of Vertices 


45 


46 


possibly at boundaries) neighbourhoods in T (0) . Such 
neighbourhoods are illustrated in Figure III.é6éd. For each 
neighbourhood, the associated vertex and its incident edges 
are dropped and the neighbourhood is re-triangularized. This 
new triangulation of T is called r(l)={ri0) ji Sel ye 2; ee nite 
where nj must be less than no since each neighbourhood now 
contains fewer triangles “Sin Jits partition (see Figure 
Tei 6 ci) 

This procedure is continued to produce t(2) from Tel) 
and) SOmone until, onlyesthe Circumscrabing triangle Ty remains. 
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FIGURE SIII46e. Triangular Subdivis1on 
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FIGURE III.6f: Triangular Subdivision TW) 
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The set of triangular subdivisions may be placed 


hierarchical tree structure as follows: 
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Analysis 
isearen! 

Et 1s? "essential to show thatthe» nerght of the above 
tree Ss tLucture™ ensuresss !ogaritimic- = "search= = comp lex itys 
Kirkpatrick shows that there exists positive constants cj 
and cz with cj,c2 > 1, such that ny] < n/c] and the number 
CE SOnsmat -each) node 1s™ bounded! by) co. |) Thus, the height =m 
OPS Chem tree sis pounded malove = bYmlOg mn, LOd cir. Rene TO mer 
at most, cj. log Nno/log cy, triangles need be examined in 
Order Com determine that: triangle 7; (0) which contains the 
POtNE 4 meeoMNCe =, 0< .C2n, swherLeuca nice constant and snaicetne 
DUNC rm Olumy CinwlCCS Slit methe m= pOlLydgOnmE>, 11cm Ol lOwSmeethat 


Kirkpatrick's search algorithm requires O(log n) operations. 


[Preprocessing] 

The preprocessing costs of Kirkpatrick's algorithm are 
as follows: Conversion of an n-vertex polygon into a 
regular graph requires O(n log. n) operation. The 
triangulation process and the hierarchical tree construction 
ate the k'th step requires, O(n, ) Operations. — Thus, thes cost 


Of the entire triangulation. and tree construction is 
m 


Mm 1 
>> O(n_) < O(n) . >, 1/ci = Ong) 
k=o 1=0 
Theretore, the total number of Operations Lequired @by “the 


preprocessing phase is O(n log n). 
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[Storage] 

ihe 'spacesireguired mofmtheatriang les methodurconsistse of 
storing n/c, tavanglées@e athe cache=tlevelifgommethne, tree 
structure. Since the height of the tree is bounded by log 
No/log cj = O(log n), requirements are bounded by 


Veg tn } 
SS No/Cis ao (11). 
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3.5 Rectangle Method 
The description of this method requires some 


preliminary delbmmitions. 


DEF 'N: A section of a polygon is defined to be a set of 


consecutive segments of the polygon. 


DEF'N: Rersection swhich 1s @monotonitCcewitherespect sLOmpoth 


BieekeandesyeCcOOrdingtesa1secal led, aeSiImplesSeCLION. 
DEF 'N=) WA maximal simple section is called a basic section. 


By = partitioning #a polygon “intow the union Cry basic 
sections (this requires determining those vertices which are 
local extremes in either the x or y coordinates), the 


processing costs of the line-crossing method (see section 
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2.5) can be significantly reduced. Let the line passing 
through a given test point and parallel to the x-axis be 
Calledmtnes test Line, For the line-crossing method, it is 
required to determine the number of intersection points of 
the test line with the boundary (with the basic sections) of 
Liew POLY OOmeandstyinge left. (oreerignt) oF the. test =poimts 
Testing to determine if the test line intersects a basic 
section, (wes shall “call)this intersection test) requares 
Simply the comparison of the y-coordinates of the section. 
Of course, for those basic sections where intersection does 
Occur, additional processing (we shall call’. thiss@the 
Orientation test) of the segments which compose the section 
is required to determine if the intersection point lies left 
One erigit Orem Lue ese point. This enhancement to the 
line=crossing/’method is the essence of the work done by 
Loomis[21], and in many applications can result in 
Significant improvements. 

ThemOrientatlone test. 1s -Simplibied by. tirst. comparing 
the x-coordinate of the point with the extreme x-coordinates 
Spy tiembasic Section, This, together with the intersection 
test, ismeduivalent to determining 1% thes tests point. lies sin 
the smallest rectangle encompassing the basic section. Only 
ifthe test sooint lies in ‘this rectangle. (called =the (section 
rectangle) does it become necessary to process the segments 


which compose the section. A geometric illustration of basic 
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section rectangles oe Ly 25 oie) tet eet OL Am DOLYGOne ise o1ven 
in Figure ITI.8a. As can be seen, the union of all the basic 
section rectangles of a polygon serves as an annulus 
approximation to the boundary of the polygon. 

The basic section rectangles can be combined pairwise 
consecutively to yield a cruder annulus approximation as 


follows: Let 


Pcees spall 


O ) 
and let 4° A SES pecta bepasnull) rectangle. Define 


(1) (on) (0) ee y 
= R2ji-1 U Beep op (akon & Sayan ay) 
| al Ge 
Sesh elas smallest Twectangle Ss whi chuencloses Ro 721) andeko aes 


Thus, re is the smallest rectangle which encloses the basic 
sections COmta lined yin Roe and Ke nese enclosing 
rectangles 1 Lome theggivenm polygon gare illustrated in 
Brgure, lt1, 6b. 

The line-crossing method can be enhanced by applying it 
FOeslia Sew BSLCUCLULC as = LoOllows. The section rectangles 
eo as ros 1 are scanned consecutively for intersections 
WLEnethemtest. line. Tie intersect 1onswith Ri GOeCSeno GMOCeiIL, 
the algorithm proceeds directly to the next section rectangle 
hinge Onlyainethe = caseuthateintersectioOnsOle ties testa ine 
with Te dOeGmEOCCUL ma Se processing meOr methenmocs1Cumpceccion 
rectangles Ho and Roa required. This essentially can 


reduce the number of rectangles to’ be processed” by 


approximately a factor of 2. 
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FIGURE III .8a: 


Basic Section Rectangles - Level 1 
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FIGURE III.8b: Section Rectangles - Level 2 


Whats  aescrrped so’ far sisitheyessence: of Burton's 5) 
rectangle method. Oneswaddittonaly construct gis given. The 
pairwise merging of rectangles iS again applied to. the 
BECLaNngGles wha to, produce theme sequence of rectangles 
ie? isles o5s=40and sovon, until’ finally a section rectangle 


k 
Rot seedectermined wilchetightily encloses the polygon. a(sec 


Bicucecm miei c ial li soma nmi ie. Oe). 
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Section Rectangles - Level 3 
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Section Rectangles - Level 4 
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FIGURE IIl1.8e: Final Enclosing Rectangle -— Level 5 
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PICUREN 111.9: Hierarchical Structure of (Section Rectangles 
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Using this hierarchical tree of the section rectangles, 


Burton applies the line-crossing method as follows: 


aay (k) 

Beginningeats the | rocteorathesehierarechys witha). elias 

: ' , k 

fevermined lie thestest polntwezatercontainedsan Rime limites 

not, then the algorithm terminates; otherwise, the section 

k-1) (k-1) 

Beccang Les Ry, uate 1) are processed, and so on down the 

hierarchy. To describe the processing of section rectangles, 

Suppose in general that it is determined at some stage of the 
ap) 

alqorithmethat. 2) °C f}) ie mele It then becomes necessary to 
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PrRocessmeachi Or sthne=sons aie bie {37 ) of 3) 
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Consider only “the section srectangie Read Since Liem plo 
ied 
cessing of £27 he identical, If the test line lies above or 


below rf2y? it cannot intersect any of the polygonal segments 


describing it. Therefore, in this case mo, further processing 
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(j-2) (j-2) ae 
processing its sons R4j-3 and R4j-2- In the remaining case 
Jel. 
when the test line intersects aera Z is not contained in 


at | 
r§i7y) (see Figure 2Til.,l0af below), then) the) parity of the 
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BLGURE VGLI.10a: Containment of Point Z in More 


than one Section Rectangle 


Of the. number of intersections of the test line with the 
polygonal segment describing et! is determined as follows. 
he Eo mec eee ome mL c™ eCndpOwnts sO heithes DOLygOndal sesegment 
describing oe Tiehmeciemeeo uli NemiltensceCLomties pOLlycgonan 
segment an odd number of times if Z lies between e] and e9 
(SCC RZ Lier GUrC mile LOG) ieeand saeey Che number son ulmesnolhier= 
WiecmCSCCu dean Ge Jenn lasso Ge et i 0a)) MeeOnCe wnempamitvesLs 


determined no. further processing of Rjj-] or its descendants 


is required. 
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Eerminates sat sthe next lower level af 216 not contained in 
Ged, G-1) . ers yee (5-1) 
co Ee a Dae om eee Lie oeCONLa INCOM nach Lic iam oni TOT 

ed. 
Dal processing continues down the respective subtree. It is 


(3-1) 


important to note that Z may be contained in both Rone iaeand 
Serer LAW I Ne PLGULe vLIT. Oa) 

The algorithm terminates early at level k-j, j<k whenever 
it is found the LPs not contained in any rectangle 
Nona eee atwthat#ilevel? The *parnrty. cf the#number of 
times the test line intersects the polygonal boundaries from 
the left then determines the solution of the point-in-polygon 
problem. Otherwise, Z is contained in one or more basic 


| °) 
sections rectangles i and additional processing of these 


rectangles is required. 


The processing of a simple section rectangles is greatly 
Simplified because the polygonal segment number which defines 
LAPS monotone with respect =toly. Binary search techniques 
Can - therefore be Used ‘to. "discriminate “quickly “the ‘point “2 
against the polygon segment. However, Burton suggests an even 
faster algorithm. The basic™Section”™ is "split into “two simple 
sections, each with an equal number of vertices (within one). 
Bacheecimplem™sectionb 1s we then tightly = enclosed ™by =a ‘simple 
section rectangle. Therefore, as in the previous hierarchical 


structure for compound sections, the basic section rectangle 
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caghtdy encloses’ sine “twomsimple ~section rectangles. The 
splitting strategy is next applied to simple sections, and so 
on until simple sections consisting of one polygonal edge 
(called primitive sections) only remains. These form the 
leaves Sol ethes total hierarchicals. tree sof compound, abasic, 
Simple and primitive section rectangles. The point-in- 
polygon algorithm previously described can be applied 
directly tomthis «total #hierarchicalmtree. Only im’ the =tesc 
point Z is contained in a primitive section, does it become 


necessary to discriminate Z against a polygonal segment. 


Analysis 

[Search] 

Let Z = /log ni.wThen total hierarchical tree contains at most 
2 ! leaves, corresponding to the primitive section 
rectangles. The total number of rectangles in the tree 
(nodes in the tree) is therefore bounded by Dig are ee 
Because processing a rectangle requires a fixed number of 
OpeGea tLOncy ald eeLiC ™ WOLSTt sCcasC miWwie emit Cm—rCeS Camm p>O 11) CuI 
contained in all the rectangles, the complexity of Burton's 
eI GOrithm ise) O(n). Indeed, Figure ill.) sprovides. jan 
illustration of a polygon where Burton's method may require 


O(n) operations to solve the point-in- polygon problem. 
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FIGURE III.11 A Worst Case Polygon (The Claw) 
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In practice, for "well-behaved" polygons, Burton's method 
PeLlOrms® signiiscantl ywbetter. The algorithm terminates at 
thatatyeveisnoiy? the ghiecnarchical @*treei ewherespnoneseoet ithe 
associated section rectangles contain the test point Z (see 
Figurespihbl sec jr Furthermore eifiiztsiscontained ingatemost 
one section rectangle, at. ‘each level, the complexity of* the 
algorithm reduces to O(log n). Indeed, because of the binary 
tree construct, it may be natural to erroneous assume this to 


be the case in general (see, for example, Fowler[10]). 


[Preprocessing] 

Thes totale numbernrofs;operationssenequaredr to,;findp-alls the 
CUENGINGEMPOLN Bsrerse O(n) Furthermore, each of O(n) section 
rectangles can be computed in fixed time. Hence, the total 


cost of preprocessing is of complexity O(n). 


[Storage] 
Since each of the O(n) section rectangles can be stored in a 
fixed number of words, the storage costs are of complexity 


O(n)s 
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3.6 The Strip Method 

The fundamental process common to these methods is the 
construction of a new polygon P* which approximates the 
Original polygon P in the following sense. Given a 
tolerance or resolution e, the distance (usually Euclidean) 
OLeanyepounceOnethe boundarysoieP eto the sboundary (Ober = sis 
bounded by e. the polygonal approximation P* is to “be 
Se leCCrec SOmeCh a alt Sw SdmpLerme ror processing spurcoses 
(Lewer pnunber Por #vertirces) |) and vinssaddition «captures =the 
chaveacter=™ sole ctiews Originals.» polygon, bP. The Smatter oF 
character capture is not discussed in depth here, because it 
has little bearing on the performance of point-in-polygon 
algorithms, except perhaps indirectly through its effect on 


the polygonal approximation P*., 


One algorithm Ele Constructing a polygonal 
aopEO ina pron) tOnwas DOLLY GON eae TnV Oe. ey Vint Tels 


given below. 


Algorithm: 
Wemekecord Viadsc san veLlLexmoL 
Set starter s = 1 
Qe wSetedestinator £°="sti 
3. Draw a line segment L joining Vg and V¢ 
Determine the maximum distance d of the points 
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Ae dia ee. ander <n 
then 
(comment: extend line segment) 
sigic if =) jest 
GOs torsten.s 
else 
(comment: line segment is complete) 
GECOLrd Veojas amgvertex. or 
GO, EO "step. 5 
Sh Abe ele 
then 
set s = t-l 
go, te, step 2 
else terminate. 

TieweeLleuty Of Sthes algorithm Ise illustrated soy eappd ying 
DesatOmcie polygon? shown eine Figures. ileal Zar The polygon 
represents the boundary of a river basin in the province of 
Alberta obtained by means of vector digitization. It was 
BLOGUCCU Mm eCOMN@e an Maps satu edamls Cal Cnn Ola) 20.70 U0 ued Cemea tamer 
resolution of approximately 50 meters. The polygon contains 
Me7omVventices » and issquites typlcalmor these polyqonsmwnich 
appear in many geographic database systems. 

The polygonal approximation P* obtained by the algorithm 
ateva tolerance ce of 250 meters: is Shown jin) Figure 111.125. 
The polygon. PF scontainseeadmstOtalweol@2iesvertices;s which jis 


approximately one percent of the original polygon P. 
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A Given Polygon P 
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PIGURE Dl w2b er OlygonalmAppLloxamacionmornes.> 
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No claims are being made that the given algorithm is 
optimal in any sense. Indeed, the algorithm described by 
Douglas and Peucker[9] is likely sto better capture the 
cartographic character of polygonal boundaries. An 
algorithm described byweeal lardi34 is somewhat more 
efficient.) And finally; theyalgorithm given by Kurozumi and 
Davis[17] yields computationally a more effective 
approximation in the sense that the polygon P* may have 
fewer vertices and yet meet the tolerance requirement. The 
given algorithm has the advantage of being simple and yet 
Maintaining some similarity with all the other algorithms. 

The polygonal approximation P* may be used to obtain an 
annulus approximation to P as follows. Move each of the 
line segments of P* outward by the distance e to obtain a 
new polygonal approximation P**, Thaseepoly gon =P eanust 
enclose the polygon P. Similarly, move each ot | the line 
segments of P* inward by the distance e to obtain a polygon 
P***Fwhich lies entirely inside P. The boundaries of P** 
and P*** describe an annulus which enclose the boundary of 
Py For the polygon P given in Figure li.i2a, the annulus 
approximation 1s illustrated in™ Figure iit l2c. In the 
figure, the annulus is partitioned into quadrangles (strnips) 


by joining pairwise the vertices of P** and DAA a 
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AScursony inspections Of Figuresiliel2c suggests that some .or 
the caricature of P is not captured by P* and consequently 
DVM oe ander. <% An enhancement in caricature capture is 
obtained by reducing the tolerance e, and thereby improving 
the accuracy of the polygonal approximation P*. Thateels, 
there is a trade-off between obtaining an approximation P* 
with improved accuracy, and increasing the number of 
vertices required to represent it. For the example under 
consideration, an improved annulus approximation of P with a 
tolerance of 50 meters (and containing 57 vertices) is given 
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FIGURE III.12d: Annulus Approximation of Polygon P 


with Tolerance of 50 Meters 
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The; questionwok tolerance jselection is a’ non-trivial 
consideration. If the application does not predicate a 
choice, One possible strategy might be to remove from P 
EHOSCEVCECICesawhich dO ™nOtecOneDlbULCeMSTUNI Ri cantly atomenrs 
definition, for example, those vertices which are nearly 
collinear, Or, tnose which Cansberateributed tor "noise". @inis 
can be achieved by obtaining successive approximation P* to 
P with ever decreasing tolerances until one is found where a 
further decrease would result in a dramatic increase in the 


number of vertices representing it. 


For faepolygoner obtained by fdigicization, ithis) strategy 
will) often lead to an approximating polygon P* which is very 
Sitelarm tO P bit which contains "Onl Vaassmnall poreron. of athe 
number of (vertices way IB For the example under 
consideration, in Figure Temes he see the annulus 
approximation of P with 50 meters and which contains only 5/7 
vertices. For a summary of experimental results for many 
ditterent. tolerances which» support™ the validity sof) “the 
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The point-in-polygon algorithm is given as follows: 

1. Apply line-crossing method to P*** 

oe ae 7a Cy Ps x 
Then terminate since ZC P 
else go to step 3 

3. Apply line-crossing method to P** 

ARGASTE Zeqe pss , 
then terminate since Z ¢ P 
else go to step 5 

oe CComment< 96*Z istcontained sin annulus) 
Apply line-crossing method to quadrangles to 
determine one (call it Q) which contains 2Z 

6. Apply line-crossing method to the line segments of P 


COntainedsine® tondetermines whethemseogmmnot. 2) Pe 


Let m be the number of vertices in P* and let Z be the 
Maximum number of line segments of P in any quadrangle. Then 
steps 1, 3 and 5 each requires O(m) operations, whereas step 
6 requirés™O(L)* Thus;eein sehes wors tecases whensrthe fonts 2 
lies in the annulus, the algorithm requires O(m) + O(2) 
operations, where mi > n. For well-behaved polygons, 
HOWeVeumeZelikely lres Outside thevannulus@and es teps soands6 
OL thesta LOOrivhnipsares not executed. The cost is therefore 


O(m) operations for the majority of points tested. 
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Note that the algorithm may be improved by combining 
Steps IM through 5." This moditrcation mowever detractc utrom 
its readability, and does not effect the asymptotic cost 
estimates given above. 

The problem of selecting an appropriate tolerance is 
solved by Ballard{3| “by means’ ofa data structure, “a Strip 
tree which provides a polygonal approximation to P for any 
bequired. sresOluelone =e. A recursive alGoritni ee, Lor 
constructing this strip tree is given. below. 

Algorithm: Strip—tree (Vpn;Vi.... V4) 

(Comment: Let (V0ypViyecss Va l= cenote™ as. polygonal 
segment) 
lis DEawsael newb thnrougnaVgjeand. Vij. 

2 Find the smallest rectangle with sides parallel to L 

WEIChE dust COVED () peaViiiessiece 7 ine 

(Comment: This) rectangle 1s called @the’ strip~cor the 

excreneee IV aalnooo oni dis 
Bis Gompute the swidth »e of the estrip was the ~cuclidean 

distance between the two sides of the rectangle which 

are parallel to L. 

4. Tie 2 

then return 

else go to step 5. 

Dis Find a vertex V, which lies on one of the two sides of 


the rectangle that are parallel to L. 
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Gr (Comment: Construct strips for segments 
WO7Viresss7 Vk) and [Vi Vi-41,2.5 «piven ) 
SUbt~ ree B(VG; Vili eee gi) 
SUDURCEN(Vi-g Vic prerers Vae 
EECUENS 
ASStlip tree for a polygon; with) Vereicesrvii,Vo7....Vn 
can be built by invoking the algorithm for the two segments 
CVV oe VON Vie, Vig iy ete VitnV ee iClLe Vineismse Lecled 
LOPbe™approximately™ "opposite GiIViPon Pe | Themroot strip 
off thesstrip cree for ‘the polygon) is definedsp to; the smallest 
rectangle» enclosing the gstrips § £00) (Vij 7V0,...) i) and 
(Vicgie sess Vries Figures smb. lob and Ile isc erliustrates the 


dara structure fora “Strips tree, of |a typical §polygonal 


Segment in Pigene Ill .lsa. 
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The polygon P can be displayed at any resolution e* by 


the following recursive algorithm described by Ballard. 


Algorithm: Display segment (T,e*) 


(Comment: Thegmeqmentewithearcotenodes 1 fis Ndisplayed 
with resolution at most e*) 

Time (T)y <ayex 

then ydisp lays stip x Ty) 

BecCuULEH 

else display segment (left son (T), e*) 

display segment’ (right son (T), e*) 


Wee Uiishie 


The point-in-polygon algorithm using the structure is 


described by Ballard as follows: 


Algoruvthm: (Vntersection) (272) 


(Comment: Computed is the parity of the number of 
times a ray from Z intersects the polygonal segments 
Within) themstGvpe.) 
Teg@vays Intersects moothmsides) Ofte ene msstri pe lethal wake 
paral le lsto #thegline beatseexstripitrees algorithm) gthen 
return (true). 
Iferay, doesenot antersect either side jofiethe: strip, \T 
thateare parallel utonlLpthenereturn (nalse)):. 
Otherwise (Comment: exclusive or) 

return [Intersection (Z, left son (T)) 
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Analysis 
[Preprocessing] 

het Tio denote they cost oryconstructing, the stwim tree 
Goma polygonal segment Mwithpeitlwvertices MIVQj,Vi,4--Vi) 
Using Ballardss algorithm. @The cost ‘oL computing themstrip 
which tightly encloses the polygonal segment, plus the cost 
Oipidentityinge@the vertex V_(0sGike aij) whicheiges onpone not 
the longitudinal sides of the strip, is bounded by ci, where 
Geis) Ja) constant ewindependant, | of 1. The strip sree 
[VO7 Vlere oe rly then becomes @avallable sas coon. as) the serip 
tree (for DOE is oes el. -teisel TIVE Aeti soo nal levee lereein 
determined. Thus the total cost of computing the strip tree 
ers NG lp aS cu ce ANS, SG) A Grey) ae. tems 

Make the severe simplifying assumptions that 1 is a 
power of2 (1. = oe say), and that k happens to be a mid- 


VErLexm thal els ko 24-1), then 
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or(2%-1l) + c2 


ohp( 20) 04 12", 


LO) 


Sincerl( 20). = 0, 16 follows sthatel (1) =3e1logeim=scOMthat,mror 
joint hemalgGOrlthteismoLmCOmplLex dt yasO(Nw LOG ml) see cmc Os te 
is claimed by Ballard to be valid in all cases. 

However, given a polygonal segment for which k = i-l, 
for this and all subsequent sub-segments (as for the segment 


in Figure III.14), 


yaa Ss 
ea7i ek a unt zngepames jo: devo raha wine na at 
147 ee th wack s380 es ate : a a ey A 


qiase sat eniseumes Te Jeon Sit ba’ X toe wae 


va0o aftg aulq .indkgse Ts, i= 
te gam 09 @oit fizpigw tf > H > vied 
con id Nd eesened BE yqiuee Satay aebie £ 
save giage . sar 4°16 saanasqebar Ang NeaNe | a 
+o gifs <B> 1eoe ee) Sldallé de selineet_asely bene 
need deed Livi. epee) Soe ives .< aaa 
s(aiuie 4¢¢ OH) Teq~es-Se demo [¢ 1408 ene euctr +a 
4 & -_itejam @ UE yhomt byhawm de 
e ai t team enobtendecs . perivivcaete sises® aftr. 
-tum ¢ af od sme peed at geilt Sane; viaee gm a) ate. 
nats |, £E=*S ima +i) oaeda)) 

's3.+ 42-Sepet's 
“See 4 res e%e > ; 
70h <4e84 oe i eclis = (257 Gag? eokleva 1B, = 1C)7 3 
gen ela in gol. oi ysdveignos 30-27 nisi wit 

,aseeo Lis ni wtiay' ad os aetia® 5 

sek = Seen, SO eee emer vty mma 
dieamieee 5u2. yo 281 aahalinenem Zi 


= 
S 
_ 


be 
iit 
a 


80 


V. 
Vi 0 
: V, 
V4 
ie _ 
—— 
ae 
a 
2% V 
iM } 
V, 
A) 
WVi-2 
M4 
Vi-3 


PLIGUREMLLL. [4 seeeoclioe Drees COnsStructron 


naLaawrdaney seer glass ‘ite on 


81 


T(i) = T(i-1l) + T(1) + ci 
= UA Gay es teleh 
=e) eed (a 1) 
Thus, LO 19> 1, 91n the worctecace se T(n) = O(n-). 


The reader should be made aware at this time, however, 
thateballavd descr besmango(n) algorithm: Lor acOnstructing a 
Strip. tree, but sronswhich «the Strips: senclose sthe;polyaonal 
segments extremely crudely. hemp Lraceicabillityeasor , the 
algorithm, therefore is severely restricted. 

[Storage] 

Since, there _,ares O(mjeastrips, thes, totals storage 
requirements are O(n). 

[Search] 

Ballard asserts that for well-behaved polygons, the 
COMPLEX bY s.Cl Mathes poOlnt—in- pol vdgon 1s ssO(OGann). Tid S 
assumes the strip tree is reasonably balanced and that one 
of the recursive calls which processes either the left or 
ELGNLMSUD—-ESEDIDS a IinesteDe oe Ole ticma | GOrlthmeex tse quick Ly 
from, tine, Gecursion. For the majority of polygons enclosed 
in geographic information systems, this assumption is indeed 
likely to be valid. However, the algorithm does require the 
Processing eOtun/2astripce tor sthespathologicaleexamples (the 
ClawWweOLVchelnemlouLree GL Wile cndethenelOLom ln etl Ca worat 


case, the algorithm is of complexity O(n). 
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See CONnc Lusion 


The following table summarizes the complexities of the 
five methods described in this chapter. 


RESULTS 
METHOD 


Preprocessing Storage Searching 


Swath O(n log2n) 


Chraan n) 
Triangle n) 
Rectangle 

S ical D 


TABLE II: Summary of Pre-Conditioned Methods 


FOG wie SUrLiciently large; the = triangles smethod a ac 
asymptotically superior to the other methods. However, it 
must be stressed that these complexities are asymptotic bounds 
only (e.g. The Swath and Triangle methods both require C log n 
operations —for searching, but. CG for =the Traangles method is 
much larger). 

FOr practical purposes, many ~Othermeconsiderations will 
effect the selection of a method. Some of these include (1) 
BiemcUmeaot tye Ote Lice datams UbUGCEUTem Or OUne Ts pULDOseSmUsIiCh 
Bower nteGcece lon Of @DOlyJOnS) aml) me UhemEL NSCS] Oa li yao leche 
method toward noise an) the polygonal data, (3) thewabulity oF 
themmethod  ctOrurepldly captures thes charactern nw Or ethe polygon, 
and (4) average costs. Further discussion of these matters 


appears in the concluding chapter. 
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CHAPTER IV 


PRACTICAL CONSIDERATIONS 


4268 Introduction 

Techniques which answer the query of whether a given 
point lies in a polygon have been discussed in the last two 
Chapters. @in practice, alregion mayticonsist@of atvarilety of 
Submegions),) }forutinstance Bolocks feotiitay town icand, it Buc 
Hequired Wwortudentit yietivatesubreg jontewhichticontainsi tartest 
DOWNIE. Two methods, rectangle enclosure and circular 
enclosure which address this problem are discussed in this 
chapter. Such methods may be considered as an extension to 
thes \pornt—-in—-polygon™ algorithm. Although the discussion 
heresy 1seginote® treatedifeinwegreat Ptdepth, ite Provides san 
OPPOLEUN Mt yamtOmgeremack FOng sa SUD) eCCLe that, does notasbelong 
properly in the previous chapters (see Aldred[1l], Nordbeck 


andakysuedt[23i)) 2 


4.2 Rectangular Enclosure 
Lew 


Xmin => WL ENC) oy es Gi KEES gly, 


Xmax Maxitxq ty XK outele © teienete y Xi) 


Yninl=] Min (yyy? eee 7 Vn) 


Ymax Max(y1];Y2reeee0eees¥n) 
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cheng th cgerectang le —Ewitheverticesmn Xia Aime) (Xin ana 
(XmaxrYmax) and (Xmax/Ymin) is the rectangle of smallest 
areagowith Jsides@yparallel kto ytheecoordinatese axis swhich 
encloses the polygon. The cost of -computing this‘ rectangle 
is, of course, O(n). Rectangles of smaller area which still 
enclosems the gepolygonts caney besnibound (see Freeman and 
chapiral ii) ) thbuteeithese iwouldslrequiresrrotations:. Such 
rotated rectangles however, are more costly to determine, 
require more storage, and adversely affect the performance 
of any subsequent processing. 

To, determine whether or notpa test point 2 = 4(x;y) 1s 
contained in a given polygon, the associated rectangle is 
first examined for containment. If 

x) SOX pre O GX OMe nae OL ey SeYninr OLY o* Ymasc 
tnengutetis fotzisand Zocannet bewscontainedtun theepolygon’: 
Otherwise, the algorithms of the previous two chapters are 


applied to resolve the point-in-polygon problem. 


423 41GrrcudariEnelosure 

A circle which encloses the polygon is determined as 
follows. First, the distances among all possible pairs of 
Verhicestofe the polygon are computed, and the largest such 
distance, d, is determined. A line segment is drawn between 
Chess two Fventicesceay 1 y anceivy which achieve this maximum 


a@istance .me wnemmidpoint exc jyo)evOnmthiis Maine sségmentascerves 


otxanfinan®! (eaten get: 
yasiteme 36 siquetaed aig mi . 
(mite 2fee weteréiaoo ofa on | [68h 
apnertes ers tarpqs, 3¢deco ett Gate 
ilj@te acter Sern sofdeme 32 snfgresace 
ans cimeerd sea) GavGet 2¢ ee pone = 
ones snaijesod Sylops: , BIpiw sgedt 2a 
e-aneieb G) yidadeas sr0n e534 as taved asilqvedoss 
s.-agrroisey, ety soeite ¢learerts Pane saperase stom % w 
ig) = 2 Sei] seo7. 4@0n- 3) 7eAIeRe saimretot a 
imetooy Lageicenae alt ..rcepieg bite ‘6. 8) 
+f ..dhsuiiedaes $o2 Benivess 


ee 


~ 


= 


Ww 


—' , ae vais ¢ 20 iecel 4 8.32. apbet 2 = 
inom ioe <tt ct Son Gjincar ed) dnaméce sere yaaa at ot 
aie eysta@ecs. owls gugilesq) any To apittlioeka CNET? %& 


wmaldorg népyinn-ai-saicg ser syvices: a3 beilg 


eobeat Ore 
=i nopi low a spears (i, aiesiv .« 
Hiaaog Ais ries soot eI) -Sriz’ <eeust 2 
Soa. J=9e4el 3249. brs » oe whos Ste tag Lg Sle to 
paeWIai nucwt ei Ihempoe, sels \beenesee as a (Sam: = 
mumrxen eit eveldie  (aviw 7 ¥ “Baa ie canta ot 5 
Pirvses snendsa dail) alA$2 ao fotteneh oe a 


if 
e 
u 


G 
» 
wi 
4 
sf 
ee 


a 
nn — 


85 


domeCiCmECeN Crem OL ICL LChomswi thee rad Miser —cd/o4 For "well- 
behavedsepolvVoOns athicmcircle mt igitlysscilrcumscntbecutne 
polygon. If it does not, the vertex V,_, outside the circle 
and furthest away from the circumference is determined. The 
distance of V,p from the centre (xc¢,yc) can then be used as 
Bad lus Ole ay newmCirclcmwithmecentremgn XorV aii LL SeecirCLe 
encompasses the polygon, but not necessarily tightly. The 
complexity of this algorithm is O(n2). 

Alternatively, a tighter circle is) obtained by passing 
a circumference through Vj, Vj and Vp. These three vertices 
can then be used as the first step in an iterative algorithm 
described briefly by Nordbeck and Rystedt[23] to find the 
circle of smallest area which encloses the polygon. 

To determine whether or not) a test point Z = (x;y) is 
Contained in given polygon, the associated) circle jwith 
Cent LemCx , Min) OOM LACTIC Tw LSM Lote ex aly LOCC gmO 1s 
containment. If (x-xc)2 + (y-yc)? > r% then it is not, and 
ZaemMus tele Outs1deisthe polygon. Otherwise, the algorithms 
Of the mprevicus (twow Chapters neeq be sapplied in order sLo0 
resolve the point-in-polygon problem. 

Circular enclosures have advantages over rectangular 
ones in that 

1. they require less storage, and 

2. they normally produce tighter enclosures for those 


polygons which are relevant to most applications. 
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On the other hand, rectangular enclosures 

1. are easier to determine (O(n) versus O(n2)), 
and 

Zee ace usually more effective ots subsequent 
computations (for example, the test for containment 
requines at (most sLour. comparisons = Versus.) three 
additions, two multiplications and one comparison 


required by a circular enclosure). 


4.4 Multiple Polygon Processing 


A common characteristic of most geographic information 
SVS Celiswels mameCapaDLilty lO ™Palulc lL Onmassdl chm Leo. Ollmmor 
interest into a number of smaller sub-regions. For the 
storage of descriptive data associated with some geographic 
point 2, it often becomes necessary to identify that 
sub-region (polygon) which contains Z. RCo Samet ees 
context that the enclosures of the two previous sections 
play a fundamental role. 

AL valgorithmss which identity sathat ssub-reqioOn es that 
contains Z require the processing of numerous sub-regions. 
Prior to invoking any particular point-in-polygon algorithm 
described in chapter Il and IIL in “order to determine 
whether or not a specific sub-region contains the point Z, 
it may be tremendously beneficial first to determine if 2 is 


contained in the encompassing rectangle or circle. If it is 
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not, then Z cannot be contained in the polygon and the much 
more costly point-in-polygon algorithm need not be 
involved. Processing then skips to the next sub-region. 

Ltt eiseeassumed@einamthe  —abover® paragraphs thatumthe 
sub-regions are listed sequentially and are processed 
accordingsto, thisslisteuntultmone Tsmeoundswiach contains =tne 
DOU Nt eZ. fo mincrease them liked 1 hocdsPorerwtinding the 
SUD-reglone= Which =“COontains 722 early in acne mvelist, a 
Parctlalordening»s ofetthe @sub=regqionsercan first” be @detined. 
For example, the sub-regions can be ordered according to the 
distances of the centres of their encompassing rectangles 
(OL pECcircles) Strom aypointezg external tom thes tocals region 
OfPpinterests The? sub-regions) closesem tostthempoint pZoyeun 
the sense implied by the ordering above, would then be 
PLOCessedPetirsit. FOme "General Bdiscussion Fotwerelated 


implementation matters the reader is referred to Aldred [1]. 
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CHAPTER V 


CONCLUDING REMARKS 


In this thesis, a variety of vector-based data 
structures for representing a polygon have been discussed 
and various algorithms for solving the point-in-polygon 
problem with respect to each of these structures have been 
presented and analysed. Analysis throughout has’ been 
confined primarily to asymptotic, worst-case analysis. The 
algorithms have been classified into three categories, 
namely, direct methods, pre-conditioned methods and multiple 


region processing methods. 


Direct methods for solving the point-in-polygon problem 
are defined to be those which operate directly on the 
polygon represented as an ordered sequence of n vertices. 
The ‘cost™oL Teach of “the direct methods ‘described i1s7 O(n) 
Pre-conditioned methods require that the polygon first be 
transformed into suitable data structures before search 
algorithms are applied. These structures permit the 
development of asymptotically faster algorithms [as fast as 
O@ogwren) | for solving the point-in-polygon problem. 
However, because of the additional preprocessing costs 
lLtypica llys O(nelogmn) lsetnmconsStLucl Inge tiescm. tLUCTULeS and 
because of the additional storage costs in representing 
them, “the” benefits of “fast search algorithms ~from a 
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practical. sense . canj. be. ~realized.. only in). a ‘static 
environment. In a dynamic environment where the boundaries 
Of thes, polygon aremesubject, toss frequent, ichanges,., but 
depending on the frequency of these changes relative to the 
frequency of searches, the use of direct methods can turn 


Out tombes superior. 


Of the direct methods, aie is found that the 
PIne=ecrosci pgm auct Od ae tOre msOlvuandeethes polnt—-in-polygon 
problem is superior in all respects. It is simpler, more 


CrLicvent pmaandmapplicablemtOuarbltranry ppOLygons. 


Of the pre-conditioned methods, the triangle method 
Will sPLeCDLOCeSS1NGisCOSESs= OfesO (ne lOguan i. SLOradeesCcosts Ob 
O(n) and search costs of O(log n) is asymptotically superior 
to the methods described. Indeed, Kirkpatrick shows that 


the method is asymptotically optimal. 


From a practical point of view, however, asymptotic 
worst-case analysis of costs is inadequate for the following 
reasons: 

Ls Trai se Juste vedaonly, fOr engsutiacientiyes large. ThUSs, 
it, is possible that the chain method whose. search costs 
are O(log2n) May be faster than the triangle method 
WHOSemSeacch BCcOStSmalemO (LOde i )mtOrma ld nmol mpnaGcd cad 


interest (say, n < 106). 
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Ze It may have no relationship to expected costs, since it 
must be valued for all cases including pathological 
ones. Thus, although the costs of the rectangle and 
Stripe methodsmare 20(in)®ifor agicertainfpolygon,init eis 
quite possiblége(weefsuspectienitisis wlikely)s thatyjthe 
average costs of these methods for random polygons and 


fon, polygonsgor mipractical/y srmoor tye syOl log nj 


To answer the questions implied in the above two points 

requires the following undertaking: 

ties The careful unbiased implementation of each of the 
algorithms. 

ae The selection of appropriate test data 
("interesting"classes and sub-classes of polygons). 

3% The selection of effective comparative evaluation 
Gr itberra. 

These matters we leave aS our primary suggestion for future 


research. 


Inpeselecring sas metnoG, shor 7am Specific application, 
factors other than the comparative performance of the point- 
in-polygon algorithms can influence the final choice. The 
data structure smuncder. Vingmecach Ol si Lhe memethOdcamay 06 
PaGticularlyeecuttablesmOl adaptable LOn msOly Ing BeOcher 


problems relevant to that application. Thus, the swath and 
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chain methods are easily adaptable to the problem of spatial 
convex inclusion and to the problem of locating a point in a 
planar subdivision induced by m straight lines. MThe chain 
and triangle methods are both suitable for solving the 
nearest neighbour (post-office) problem, which is useful in 
awe VaLLlecy = Of mapplicatronse »incivairng “surve Plance= eand 
Gracin. In addition, the triangle method is appropriate 
in applications involving triangular finite element models 
(eg. some terrain models). The rectangle and strip methods, 
are conducive to solving problems of spacial analysis, 
including the union and intersection of polygons. 
Furthermore, because these last two methods yield 
approximations to the boundaries of a polygon, they can be 
used to capture the caricature of polygons at various 
resolutions. This is useful in numerous applications 


including graphic ‘display ‘systems. 
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